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I analyze common problems of interlayer tunneling in Bi-2212 mesa structures, such as self-heating and 
nonuniformity of junctions. Numerical simulations have shown that self-heating does not mask the temperature 
dependence of the superconducting gap. Major problems can be avoided by decreasing mesa sizes. 
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1. Introduction 

Layered high-T c superconductors (HTSC) rep- 
resent natural stacks of atomic scale intrinsic 
Josephson junctions (IJJ's). Indeed, interlayer 
spacing of 15.5 A was estimated from a periodic 
Fraunhofcr modulation of Fiske steps in Bi2212 
mesas |5J]. IJJ's with their record-large I c R n val- 
ues ~ lOmV are attractive candidates for cry- 
oelectronics. Interlayer tunneling spectroscopy 
has become a powerful tool for fundamental stud- 
ies of HTSC However, interlayer tunneling 
has several problems, such as defects, nonunifor- 
mity of IJJ's in the mesa and self- heating or non- 
equilibrium quasiparticle injection during trans- 
port measurements. 

Here I analyze typical consequences of those 
problems. Numerical simulations have shown 
that self-heating does not mask temperature de- 
pendence of the superconducting gap close to T c . 
Major problems could be effectively avoided by 
decreasing the in-plane size of mesas due to a pro- 
portional decrease of defects and self-heating. 

Fig. 1 shows experimental current-voltage 
characteristics (IVC's), normalized per IJJ, 
(three probe configuration) for an intercalated || 
HgBr 2 -Bi2212 mesa 10 x 20 fim 2 at 3.9K < T < 
80K > T c ~ 77K. IVC's reveal a behavior typ- 
ical for SIS tunnel junctions. The knee in the 
IVC's represents the sum-gap voltage V s = 2A/e, 
which vanishes with approaching T c pj and is sup- 



pressed by magnetic field j3J. Temperature de- 
pendence of the superconducting gap, A, is shown 
in inset together with the BCS curve. 

2. Self-heating 

An analytic solution of the self-heating problem 
for a circular mesa yields fl: 
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where T is the base temperature, a is the ra- 
dius of the mesa, q is the dissipated power den- 
sity and n a b is the in-plane thermal conductivity. 
For Bi2212, self-heating may become large at low 
T due to a small K a b- From Fig. 1 it is seen 
that there is back-bending of IVC's at the sum- 
gap voltage at low T, due to self-heating and/or 
nonequilibrium quasiparticle injection. 

Fig. 2 shows results of a numerical simula- 
tion of self-heating in a circular mesa a — 5/im, 
containing TV — 10 IJJ's with parameters typi- 
cal for optimally doped Bi2212 |g{ (p c = 30ficm, 
A(T = 0) = 30 meV, K ab (T) ~ 10~ 3 T (K) W 
cm - 1 K _1 , T c = 95 K). Simulations were made 
in the following manner: First, a set of IVC's(T) 
without self-heating was calculated, see Fig. 2 
a), using standard SIS tunneling equations and 
a model A(T) dependence, shown by the solid 
line in Fig. 2 d). Calculations were made for d- 
wave superconductors with coherent, directional 
H tunneling, providing a qualitative resemblance 
to experimental data in Fig. 1. Next, the current 
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Figure 1. Experimental IVC's of a HgBr 2 -Bi2212 
mesa at 3.9K < T < 80K. Inset shows the tem- 
perature dependence of the superconducting gap. 



was swept for a fixed T and equilibrium V and T 
were obtained by numerically solving Eq.(l). Re- 
sulting IVC's for several T are shown in Fig. 2 b) . 
It is seen that self-heating causes back-bending of 
the IVC's, due to an increase of the effective mesa 
temperature, as shown in Fig. 2 c). 

The overheating is strongest at low T both be- 
cause of small k and large A. The "measured" 
A(T), obtained from the maximum dl/dV for in- 
creasing current, is shown in Fig. 2d). It is 
smaller than the virgin A(T), however retains the 
general shape. Indeed, at low T, despite large 
self-heating, the variation of A is small due to 
a flat A(T) dependence, while close to T c , self- 
heating vanishes as q oc A 2 — > 0. Similar ar- 
guments apply for nonequilibrium quasiparticle 
injection (/ oc A — > at T c ). 

3. Non-uniform junctions 

Fig. 3 a) shows quasiparticle branches in IVC's, 
(four probe configuration), for a Bi2212 mesa at 
T = 4.2K. Each branch represents switching of 
an additional IJJ into a resistive state. If all IJJ's 
in the mesa were identical, the multi-branch pat- 
tern would have perfect periodicity g], while in 
Fig. 3 a) a "ghost" branch "s" appears leading 
to a doubling of the branch structure. The to- 
tal set of branches corresponds to four identical 
"w" junctions and one "s" junction. Junction "s" 
has a smaller V(I) and may either have smaller A 
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Figure 2. Numerical simulation of self- heating in 
a Bi2212 mesa: a) and b) IVC's without and with 
self- heating at different Tq; c) Temperature of the 
mesa along the IVC's from b); d) Virgin (solid 
line) and measured (symbols) A(T). 



(the IVC is squeezed along V— axis) or have larger 
critical current, I c (the IVC is stretched along 
/—axis). The latter assumption provides excel- 
lent fit to the whole branch set if I c (s) ~ 2.1I c (w), 
as shown by dashed lines in Fig. 3 a) for Is, Iw 
and ls+lw branches. 

A decisive argument is provided from the anal- 
ysis of sub-gap phonon structures Q in IVC's, 
shown in Fig. 3 b). For the lw branch, three 
well recognizable phonon structures at ~ 6.3, 8.1, 
and 12 mV, are marked in Fig. 3 b). Those fea- 
tures are repeated at exactly the same current for 
the 2w branch, indicating excellent uniformity of 
w-junctions. However, when the "s" branch is 
involved, three additional maxima in dV/dl ap- 
pear. Those features correspond to exactly the 
same phonon energies but occur at ~ 2.1 times 
larger currents. Therefore, junction "s" has 2.1 
times larger I c but the same A as "w". Most 
probably this is an underetched bottom junction 
with a larger area (note that in a four probe con- 
figuration the top junction is not at the surface 
of the crystal). 

Fig. 4 a) shows experimental dI/dV(V) curves 
(three probe configuration) for two Bi2212 mesas 
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Figure 3. a) Quasiparticle IVC's of a Bi2212 mesa 
with a " ghost" branch " s" . Dashed lines repre- 
sent a fit assuming that "s" has 2.1 times larger 
area than " w" . b) Sub-gap phonon structures for 
different branches. In "s" the same phonon struc- 
tures appear at ~ 2.1 times larger current. 
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Figure 4. a) Experimental dI/dV(V) curves for 
two Bi2212 mesas from the same crystal, b) 
Simulated IVC's for a single junction (solid line) 
and for a mesa with three slightly different junc- 
tions (symbols). Corresponding dl/dV curves are 
shown in the inset. 



on the same single crystal. Both mesas contain 
N = 7 IJJ's. It is seen that one mesa shows a 
single sharp peak in conductance, while the other 
reveals several maxima with smaller amplitudes. 

Fig. 4 b) presents a numerical simulation of 
IVC's of a mesa with nonuniform IJJ's. The solid 
line represents the IVC for a d-wave supercon- 
ductor with coherent tunneling (without direc- 
tionality cf. with Fig. la)). The corresponding 
tunneling conductance curve has a sharp maxi- 
mum at the maximum sum-gap voltage (in con- 
trast to usually considered incoherent d-wave tun- 
neling). Symbols represent the IVC per junction 
for a mesa with N = 3, where the bottom IJJ 
has 20% larger area and I c , the surface IJJ has 
20% smaller A and I c than the middle IJJ (solid 
line). The total IVC is strongly smeared and the 
corresponding dI/dV(V) curve shows three small 
maxima, in qualitative agreement with Fig. 4. a. 

4. Conclusions 

I have studied the influence of self-heating and 
non-uniformity of junctions for inter layer tunnel- 
ing in HTSC mesas. Numerical analysis shows 
that interlayer tunneling provides reliable infor- 
mation about the superconducting gap even in the 
presence of self-heating. This is important in view 
of the highly debated issue of whether A(T) — * 
at T c [|| or not 0. Apparently, self- heating can 



not explain the observed strong drop of A(T — * 
T c ), if A were completely T— independent. 

Non-uniformity of IJJ's deteriorates interlayer 
tunneling characteristics. The only way of dis- 
tinguishing genuine and defect-induced properties 
is to check the reproducibility. E.g. in case of 
Fig. 4, similar single peak curves were observed 
for ~ 80% of mesas, while the multi-peak struc- 
ture (when present) was different for each mesa. 
By decreasing the in-plane size and the number 
of junctions in the mesa, major problems can be 
effectively avoided due to a proportional decrease 
of defects and self-heating. 
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